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Abstract. A description is given of BTeV, a proposed program at the 
Fermilab collider sited at the CO intersection region. The main goals are 
measurement of mixing, CP violation and rare decays in both the b and 
charm systems. The detector is a two- arm-forward spectrometer capable 
of triggering on detached vertices and dileptons, and possessing excellent 
particle identification, electron, photon and muon detection. 

I INTRODUCTION 

BTeV is a Fermilab collider program whose main goals are to measure mixing, 
CP violation and rare decays in the b and c systems. Using the new Main injector, 
now under construction, the collider will produce on the order of 4 x 10 11 b hadrons 
in 10 7 sec. of running. This compares favorably with e + e~ colliders operating at the 
T(4S) resonance. These machines, at their design luminosities of 3 x 10 33 cm _2 s _1 
will produce 6 x 10 7 B mesons in 10 7 seconds [1]. 

II IMPORTANCE OF HEAVY QUARK DECAYS 

The physical point-like states of nature that have both strong and electroweak 
interactions, the quarks, are mixtures of base states described by the Cabibbo- 
Kobayashi-Maskawa matrix: [2] 

d'\ ( V ud V us V ub \ ( d \ 

s' \ = \ V cd V cs V cb s (1) 
V J \ Vu V ts V tb ) \b ) 

The unprimed states are the mass eigenstates, while the primed states denote the 
weak eigenstates. A similar matrix describing neutrino mixing is possible if the 
neutrinos are not massless. 

There are nine complex CKM elements. These 18 numbers can be reduced to four 
independent quantities by applying unitarity constraints and using the fact that the 
phases of the quark wave functions are arbitrary. These four remaining numbers are 
fundamental constants of nature that need to be determined experimentally, like 
any other fundamental constant such as a or G. In the Wolfenstein approximation 
the matrix is written as [3] 
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( l-A 2 /2 A AX 3 (p-i v ) 

-A 1 - A 2 /2 AX 2 

\ AX 3 (1 -p-irj) -AX 2 1 



(2) 



The constants A and A have been measured [4]. 

The phase 77 allows for CP violation. CP violation thus far has been seen only in 
the neutral kaon system. If we can find CP violation in the B system we could see 
if the CKM model works or perhaps go beyond the model. Speculation has it that 
CP violation is responsible for the baryon-antibaryon asymmetry in our section of 
the Universe. If so, to understand the mechanism of CP violation is critical in our 
conjectures of why we exist [5]. 

Unitarity of the CKM matrix leads to the constraint triangle shown in Fig. 1. 
The left side can be measured using charmless semileptonic b decays, while the 
right side can be measured by using the ratio of B s to Bd mixing. The angles can 
be found by measuring CP violating asymmetries in hadronic B decays. 
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FIGURE 1. The unitarity triangle shown in the p — r\ plane. The left side is determined by 
measurements of b — ► u/b — ► c and the right side can be determined using mixing measurements 
in the B s and Bd systems. The angles can be found by making measurements of CP violating 
asymmetries in hadronic B decays. 

The current status of constraints on p and 77 is shown in Fig. 2. One constraint 
on p and r\ is given by the K° L CP violation measurement (e) [6] , where the largest 
error arises from theoretical uncertainty. Other constraints come from current 
measurements on V u b/V c b, and Bd mixing [4]. The widths of both of these bands 
are dominated by theoretical errors. Note that the errors used are ±lcr. This shows 
that the data are consistent with the standard model but do not pin down p and 

7). 

It is crucial to check if measurements of the sides and angles are consistent, i.e., 
whether or not they actually form a triangle. The standard model is incomplete. 
It has many parameters including the four CKM numbers, six quark masses, gauge 
boson masses and coupling constants. Perhaps measurements of the angles and 
sides of the unitarity triangle will bring us beyond the standard model and show 
us how these parameters are related, or what is missing. 
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FIGURE 2. The regions in p — rj space (shaded) consistent with measurements of CP violation 
in K1 decay (e), V u b/V c b in semileptonic B decay, B° d mixing, and the excluded region from limits 
on B° mixing. The allowed region is defined by the overlap of the 3 permitted areas, and is where 
the apex of the CKM triangle sits. The bands represent ±lcr errors. The error on the Bd mixing 
band is dominated by the parameter fg. Here the range is taken as 240 > fg > 160 MeV. 

Ill THE MAIN PHYSICS GOALS OF BTEV 



A Physics Goals For B's 



Here we briefly list the main physics goals of BTeV for studies of the b quark. 

• Precision measurements of B s mixing, both the time evolution, x s , and the 
lifetime difference, Ar, between the positive CP and negative CP final states. 

• Measurement of the "CP violating" angles a and 7. We will use B° — > 7r + 7r~ 
for a [7] and measure 7 using several different methods including measuring the 
time dependent asymmetry in B° — > DfK^ 1 , and measuring the decay rates B + — > 
D°K + and B~ -> D K + , where the D° can decay directly or via a doubly Cabibbo 
suppressed decay mode [8,9]. 

• Search for rare final states such as K[x + ^ and irfi + fi~ which could result from 
new high mass particles coupling to b quarks. 

• We assume that the CP violating angle (3 will have already been measured by 
using B° — > i/)K s , but we will be able to significantly reduce the error. 
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B The Main Physics Goals for charm 



According to the standard model, charm mixing and CP violating effects should 
be "small." Thus charm provides an excellent place for non-standard model effects 
to appear. Specific goals are listed below. 

• Search for mixing in D° decay, by looking for both the rate of wrong sign 
decay, rp, and the width difference between positive CP and negative CP eigenstate 
decays, Ar. The current upper limit on ro is 3.7 x 10 -3 , while the standard model 
expectation is r D < 1CT 7 [10]. 

• Search for CP violation in D°. Here we have the advantage over b decays 
that there is a large D* + signal which tags the initial flavor of the D° through the 
decay D* + — > n + D°. Similarly D*~ decays tag the flavor of initial D . The current 
experimental upper limits on CP violating asymmetries are on the order of 10%, 
while the standard model prediction is about 0.1% [11]. 

• Search for direct CP violation in charm using D + and Df decays. 

• Search for rare decays of charm, which if found would signal new physics. 

C Other b and charm Physics Goals 

There are many other physics topics that can be addressed by BTeV. A short 
list is given here. 

• Measurement of the bb production cross section and correlations between the 
b and the b in the forward direction. 

• Measurement of the B c production cross section and decays. 

• The spectroscopy of b baryons. 

• Precision measurement of V c b using the usual mesonic decay modes and the 
baryonic decay mode A& — > A c £~v to check the form factor shape predictions. 

• Precision measurement of V u b/V c b using the usual mesonic decay modes. 

• Measurement of the cc production cross section and correlations between the 
c and the c in the forward direction. 

• Precision measurement of V c d and the form factors in the decays D — > nl + is 
and D — > pl + v. 

• Precision measurement of V cs and the form factors in the decay D — > K*l + v. 

IV CHARACTERISTICS OF HADRONIC B 

PRODUCTION 

It is often customary to characterize heavy quark production in hadron collisions 
with the two variables p t and rj. The latter variable was invented by those who 
studied high energy cosmic rays and is assigned the value 



r) = -Zn (tan (0/2)) , 
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(3) 



where 9 is the angle of the particle with respect to the beam direction. 

According to QCD-based calculations of 6-quark production, the 6's are produced 
"uniformly" in i] and have a truncated transverse momentum, p t , spectrum, charac- 
terized by a mean value approximately equal to the B mass [12]. The distribution 
in i] is shown in Fig. 3. 



B hadrons at the Tevatron 




1i T) 

FIGURE 3. The B yield versus ij (left). (3j of the B versus r\ (right). Both plots are for the 
Tevatron. 



There is a strong correlation between the B momentum and r/. Shown also in 
Fig. 3 is the f3^y of the B hadron versus 77. It can clearly be seen that near rj of 
zero, /?7 1/2, while at larger values of \rj\, f3^ can easily reach values of 6. This 
is important because the observed decay length increases with (3j and furthermore 
the absolute momenta of the decay products are larger allowing for a suppression 
of the multiple scattering error. 

The "flat" 77 distribution hides an important correlation of bb production at 
hadron colliders. In Fig. 4 the production angle of the hadron containing a b quark 
is plotted versus the production angle of the hadron containing a b quark according 
to the Pythia generator. There is a very strong correlation in the forward direction 
(the direction of the p beam at 0°-0°), where both B and B hadrons are going 
in the same direction. The same strong correlation is present in the p direction. 
This correlation is not present in the central region (near 90°). By instrumenting 
a relatively small region of angular phase space, a large number of bb pairs can be 
detected. Furthermore the _B's populating the two "forward" regions have large 
values of /3j. 

Charm production is similar to b production but more copius. Current theoretical 
estimates are that charm is 1-2% of the total pp cross section. 

Table 1 gives the relevant Tevatron parameters. We expect to eventually run 
at a luminosity of 2 x 10 32 cm~ 2 s _1 . A machine design that holds the luminosity 
constant at this value, called "luminosity leveling," has been developed. We plan 
to adopt this design. 
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FIGURE 4. The production angle (in degrees) for the hadron containing a b quark plotted 
versus the production angle for the hadron containing a b quark. 

V THE EXPERIMENTAL TECHNIQUE: A FORWARD 

TWO-ARM SPECTROMETER 

A sketch of the apparatus is shown in Fig. 5. The two- arm spectrometer fits in 
the expanded CO interaction region, which is being excavated. The magnet that 
we will use, called SM3, exists at Fermilab. The other important parts of the 
experiment include the vertex detector, the RICH detectors, the EM calorimeters 
and the muon system. 

The angle subtended is approximately ±300 mr in both plan and elevation views. 
The vertex detector is a multi-plane pixel device that sits inside the beam pipe. 
The baseline design has 31 stations with triplets in each station. The detector is 

TABLE 1. The Tevatron as a b and c source 
for CO in Run II. 



Luminosity 

b cross section 100 /ib 

# of 6's per 10 7 sec 4 x 10 11 

b fraction 0.2% 

c cross section > 500 /ib 

Bunch spacing 132 ns 

Luminous region length a z — 30 cm 
Luminous region length a x a y = ss 50 fiin 

Interactions/crossing < 2 > 
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BTeV Layout 



Pixel Detector 
31 Stations of (yxy) 

Dipole Magnet Steel 
Field Integral '. 2.6 T—m 
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Beam Pipe 



Wire Chambers 
Aperature ±300 mr 
RICH: flourine gas radiator, 
areogel preradiator 

EM calorimeter 

Muon absorber and Toroid 



FIGURE 5. Sketch of the BTeV spectrometer. 



12 m 































n 


















\ 


— 

1 










= 









sketched in Fig. 6. Our new baseline detector has a square hole, 12 mm x 12 mm 
around the beam, instead of a 12 mm gap between top and bottom halves. (Some 
of our simulations have been done with the detector with the gap, called "EOI," 
and some have been done with the "square hole.") The triggering concept is to 
pipeline the data and to detect detached b or c vertices in the first trigger level. 
The vertex detector is put in the magnetic field in order to insure that the tracks 
considered for vertex based triggers do not have large multiple scattering because 
they are low momentum. 



The RICH detector [13] has a gas radiator, either C4F10 or C5F12, and mirrors 
that focus the Cherenkov light onto photodectors situated outside of the fiducial 
volume of the detector. This system will provide K/tt separation in the momentum 
range between 3 and 70 GeV/c. To resolve protons from kaons below the kaon 
threshold of 9 GeV/c, a thin aerogel radiator may be placed in front of the gas 
volume. The same photon detector would be utilized. 



The muon detector consists of position-measuring chambers placed around and 
between an iron slab followed by another slab used as a magnetized toroid. This 
system is used both to trigger on final states with dimuons and to identify muons 
in the final analysis. 
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Details of Reference Pixel Detector 
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FIGURE 6. Layout of the BTeV pixel detector. There are 31 stations of triplets with the 
narrow pixel dimension in the bend plane. The most recent version replaces the 12 mm gap 
between detector planes with a 12 mm x 12 mm square hole centered on the beam. 



VI SIMULATIONS 

A Introduction 

We have developed several fast simulation packages to verify the basic BTeV 
concepts and aid in the final design. The trigger simulations, discussed below, 
are done with full pattern recognition. The input consists only of hits which are 
smeared by their resolution. To simulate backgrounds in the final physics analysis, 
we use a fast simulation which simulates track resolutions but not the pattern 
recognition. This is done because we have to simulate backgrounds in processes with 
branching ratios in the 10 -5 -10~ 6 range and we cannot afford the computer time. 
The key program in our system is MCFast [14]. Charged tracks are generated and 
traced through different material volumes including detector resolution, multiple 
scattering and efficiency. This allows us to measure acceptances and resolutions in 
a fast reliable manner. 



B Trigger Simulations 

We simulate the trigger using the baseline pixel detector shown in Fig. 6. The 
triplet stations each provide a three-dimensional space point as well as a track 
direction mini-vector. This is useful for fast pattern recognition. The trigger simu- 
lations are carried out by doing the complete pattern recognition from the hits left 
in the detector by tracks and converted photons. 

Our baseline trigger algorithm works by first determining the main event vertex 
and then finding how many tracks miss this vertex by na, where a refers to the 
impact parameter divided by its error. Furthermore, a requirement is then placed 
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on the track momentum in the bend plane, p y , as determined on line. The pre- 
liminary results of simulating this algorithm are shown in Fig. 7 (right) for a cut 
p y > 0.5 GeV/ c [15] . The choice of the number of tracks and the impact parameter 
requirement must eventually be optimized, but what is shown here is the efficiency 
for accepting light quark events (u, d, and s) for various choices on the number 
of tracks (curves) and the size of their required impact parameter divided by the 
error in impact parameter. The efficiency for accepting B° — > 7r + 7r~ is shown in 
the left side. Here the efficiency is given after requiring that both tracks are in the 
spectrometer and accepted for further analysis. For a "typical" no cut of 3 and 
track requirement of 2, the ti + tc~ trigger efficiency is about 45%, while the light 
quark background has an efficiency of about 0.8%. Note, that we do not consider 
c to be a background in this experiment. For a "typical" charm reaction the same 
trigger gives substantially less than 1% efficiency on charm, while the efficiency for 
two-body charm decays is approximately 1%. 



4 tracks " 



Impact parameter cut in a 



Impact parameter cut in c 



FIGURE 7. (left) Trigger efficiency for B° — > it + it~ for pion tracks in the spectrometer, (right) 
Trigger efficiency for light quark events. The ordinate gives the choice of cut value on the impact 
parameter in terms of number of standard deviations (a) of the track from the primary vertex. 
The curves show the effect of requiring different numbers of tracks. 



The trigger efficiency on states of interest is correlated with the analysis criteria 
used to reject background. These criteria generally are focussed upon insuring that 
the B decay track candidates come from a detached vertex, that the momentum 
vector from the B point back to the primary interaction vertex, and that there 
are no other tracks consistent with the B vertex. When the analysis criteria are 
applied first and the trigger efficiency evaluated after, the trigger efficiency defined 
in this manner is larger. In Fig. 8 we show the efficiency to trigger on B s — > if)K*°, 
ip — > K*° — > K~7T + using the tracking trigger only for events with the 

four tracks in the geometric acceptance, and the efficiency evaluated after all the 
analysis cuts have been applied. Here the trigger efficiencies for 3a and 2 tracks are 
67% for events with all 4 tracks in the geometrical acceptance and 84% on events 
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after all the analysis cuts have been applied. 




FIGURE 8. Trigger efficiency of B s -> tpK*°, ip -> fi+ >~ , A'*° -> A'-tt+ for tracks in the 
geometric Acceptance (left) and after all analysis cuts (right). The abscissa gives the value of the 
impact parameter in terms of number of standard deviations (a) of the track from the primary 
vertex. The curves show the effect of requiring different numbers of tracks. All tracks are required 
to have at least 0.5 GeV/c momentum in the bend plane. 

At the BTeV design luminosity of 2xl0 32 cm~ 2 s~ 1 there is an average of two in- 
teractions per beam crossing. The interactions are spread out over the long (a=30 
cm) interaction region. The trigger must not fire merely due to the presence of two 
nearby interactions. To insure this we have imposed a requirement that the max- 
imum impact parameter of a track not be larger than 2 mm. The yield for events 
containing &B° —> 7r + 7r~ decay as a function of luminosity is shown in Fig. 9 (left). 
Here we do not want the trigger rate to increase as a function of luminosity, even 
though this means that the efficiency on this rare final state increases. Therefore, 
a linear rise would be ideal. On the right side we show the probability to trigger 
on light quark background. We would like this to remain constant with increasing 
luminosity. No increase occurs up to a luminosity of 10 32 cm -2 s -1 , after which the 
probability for this particular trigger condition increases mildly. However, the first 
level trigger rate is clearly much lower than the 1% we require until we exceed a 
luminosity of ~ 3x 10 32 cm _2 s~ 1 . 



C Measurement of the CP violating asymmetry in 

B° -> 7T+7T- 

The trigger efficiency for this mode has already been discussed. For the B° — ► 
7r + 7r~ channel BTeV has compared the offline fully reconstructed decay length 
distributions in the forward geometry with that of a detector configured to work in 
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FIGURE 9. Luminosity dependent trigger efficiencies for B° — > ir + n (left), and light quark 
events (right). The abscissa is in units of 10 32 cm _2 s _1 . 

the central region. The left side plot in Fig. 10 shows the B momentum distribution 
and decay distance error as a function of b momentum. 

The right plot in Fig. 10 shows the normalized decay length expressed in terms 
of L/a where L is the decay length and a is the error on L for the B° — ► 7r + 7r~ 
decay [16]. 

The forward detector clearly has a much more favorable L/a distribution, which 
is due to the excellent proper time resolution. The ability to keep high efficiency 
in the trigger and analysis levels and devastate the backgrounds relies primarily on 
having the excellent L/a distribution shown for the forward detector. 

For this analysis L/a is required to be > 15. Each pion track is required to miss 
the primary vertex by a distance /error > 5a and the B° candidate is required to 
point back to the primary with a distance /error < 2a. Furthermore, each track 
must be identified as a pion and not a kaon in the RICH detector. Without particle 
identification it is impossible to distinguish B° — > 7r + 7r~ from the combination of 
B° -> KV, B s -> K+K~ and B s -> KV, as is shown on Fig. 11. Here 
B(B° -> KV) is taken as 1.5 x 10" 5 and B(B° -> n + n-) is taken as 0.75 x 10~ 5 , 
from recent CLEO measurements [17]. The B s decay into K + K~ is assumed to 
have the same rate as the B° decay into K n^, and the B s decay into K ir^ is 
assumed to have the same rate as the B° decay into it + it~ . 

Using the good particle identification, BTeV predicts that they can measure the 
CP violating asymmetry in 7r + 7r~ to ±0.013 as detailed in Table 2. 

D Flavor tagging 

We have assumed a flavor tagging efficiency of 10%. Actually our studies show 
that we probably can achieve a higher efficiency. The usual definitions are: N is 
the number of reconstructed signal events, Nr is the number of right sign flavor 
tags, Nw is the number of wrong sign flavor tags, e is the efficiency (given by 
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FIGURE 10. (left) The B momentum distribution for events in the detector acceptance (dashed 
line) and the error on the spatial distance of the tt + tt~ decay vertex from the primary vertex (solid 
points), (right) Comparison of the L/a distributions for the decay B° — > 7r + 7r~ in central and 
forward detectors produced at a hadron collider with a center of mass energy of 1.8 TeV. 



[N R + N w ]/N) and D is the dilution (given by [N R - N W ]/[N R + N w )). The 
quantity of interest is eD 2 which when multiplied by N gives the effective number 
of events useful for the calculation of an asymmetry error. 

We have investigated the feasibility of tagging kaons using a gas Ring Imaging 
Cherenkov Counter (RICH) in a forward geometry and compared it with what is 
possible in a central geometry using Time-of-Flight counters with good, 100 ps, 
resolution. For the forward detector the momentum coverage required is between 
3 and 70 GeV/c. The lower momentum value is determined by our desire to tag 
charged kaons for mixing and CP violation measurements, while the upper limit 
comes from distinguishing the final states 7r + 7r~ , K + n~ and K + K~ . The momen- 
tum range is much lower in the central detector but does have a long tail out to 
about 5 GeV/c. Either C4F10 or C5F12 have pion thresholds of about 2.5 GeV/c. 
The kaon and proton thresholds for the first gas are 9 and 17 GeV/c, respectively. 

The BTeV RICH was simulated using the current CO geometry with MCFast. 
Fig. 12 shows the number of identified kaons plotted versus their impact parameter 
divided by the error in the impact parameter for both right-sign and wrong-sign 
kaons. A right-sign kaon is a kaon that properly tags the flavor of the other B 
at production. We expect some wrong-sign kaons from mixing and charm decays. 
Many others just come from the primary. A cut on the impact-parameter standard- 
deviation plot at 3.5cr gives an overall eD 2 of 6%. This number is reduced to 5% 
because of bb mixing [19]. Without the aerogel preradiator to distinguish protons 
from kaons below threshold we would experience an additional reduction down to 
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FIGURE 11. Invariant mass distributions of all B — > h + h~ final states, where /i denotes either 
a pion or kaon, and the mass is computed assuming that both tracks are pions. The plot on 
the left shows the individual background channels and the one on the right shows the sum of all 
channels properly normalized (see text) to the ir + ir~ signal. 



about 4%. These numbers are for a perfect RICH system. Putting in a fake rate 
of several percent, however, does not significantly change the conclusion. 

The simulation of the central detector gives much poorer numbers. In Fig. 12 
eD 2 for both the forward and central detectors are shown as a function of the kaon 
impact parameter (protons have been ignored). It is difficult to get eD 2 of more 
than 1.5% in the central detector. 

Now let us consider other tags. We have simulated muon and electron flavor tags 
in our system. Although this technique is very useful at e + e~ colliders operating at 
the T (4S) , it is less useful here because it is difficult to distinguish leptons from the 
b — > c — > £ + decay from the primary leptons from the b quark decay. Our estimates 
are given in Table 3 along with those for a central detector. 

The two other methods considered are "jet charge" and "same side" tagging 
(sst). We have not yet studied sst, which is using the charge of a track closest in 
phase space to the reconstructed B. However, CDF has measured eD 2 for it to 
be (1.5±0.4)% and take 2% as their future projection using an improved vertex 
detector. We have studied jet charge, which involves taking a weighted measure 
of the charge of the tagging b jet. However, we incorporate information on the 
detachment of the tracks to help us define the jet. CDF extrapolates to 3% while 
we expect 6.5%. Table 3 summarizes our projected tagging efficiencies. 
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TABLE 2. Numbers entering into the accuracy in measuring the CP violating 
asymmetry in B° — > 7r + 7r~. 



Quantity 


Value 


Cross section 


100 fib 


Luminosity 


2 x 10 32 


# of S°/2xl0 7 s, C leveled 


2.8 x 10 11 




0.75 x 10~ 5 


Reconstruction efficiency 


0.08 


Triggering efficiency (after all other cuts) 


0.72 


# Of 7T + 7T~ 


128,000 


eD 2 for flavor tags (K ± , same + opposite sign jet tags) 


0.1 


# of tagged 7r + 7r~ 


12,800 


Signal/Background 


0.9 


Error in asymmetry (including background) 


±0.013 



TABLE 3. The projected flavor tagging efficiencies for a 
central detector similar to CDF and BTeV in units of eD 2 . 





A* 




e± 


SST 


Jet Charge 


Sum 


BTeV 


5% 


1.6% 


1.0% 


>2% 


6.5% 


>10% 


Central 


0% 


1.0% 


0.7% 


2% 


3% 


«5% 



E Measurement of B s mixing 

BTeV has studied the feasibility of measuring the B s mixing parameter x s = 
Am s /T s . This measurement is key to obtaining the right side of the unitarity 
triangle shown in Fig. 1. Current limits on B s mixing from LEP give x s > 15 [20]. 
Recall that for Bd mesons, x = 0.73. The oscillation length for B s mixing is at 
least a factor of 20 shorter and may approach a factor of 100! 

BTeV has investigated two final states that can be used. The first, tpK* , ip — > 
ji + ji~ and K*° — > K~n + , has several advantages. It can be selected using either 
a dilepton or detached vertex trigger. Backgrounds can be reduced in the analysis 
by requiring consistency with the ip and K*° masses. Furthermore, it should have 
excellent time resolution as there are four tracks coming directly from the B decay 
vertex. The resolution in proper time is 42 fs. The one disadvantage is that the 
decay is Cabibbo suppressed, the Cabibbo allowed channel being ip(f> which is useless 
for mixing studies. The branching ratio therefore is predicted to have the low value 
of 8.5 x KT 5 . 

The time distributions of the unmixed and mixed decays are shown in Fig. 13, 
along with a calculation of the likelihood of there being an oscillation as determined 
by fits to the time distributions. Background and wrong tags are included. The 
fitting procedure correctly finds the input value of x s = 40. The danger is that a 
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FIGURE 12. (left) L/a distributions in BTcV for K ± impact parameters for right sign (un- 
shaded) and wrong sign (shaded) tags, (right) Overall eD 2 values from kaon tagging for a forward 
detector containing a flourine based RICH versus a central detector with 100 ps time of flight 
resolution as a function of kaon impact parameter in units of L/a. (Protons and bb mixing have 
been ignored in both cases.) 



wrong solution will be found. The dashed line shows the change in likelihood cor- 
responding to 5 standard deviations. If our criterion is that the next best solution 
be greater than 5a, then this is the best that can be done with one year's worth of 
data in this mode. Once a clean solution is found, the error on x s is quite small, 
being ±0.15 in this case. 

BTeV has also investigated the D+n~ decay of the B s , with Df — > 0tt + . It turns 
out that the lifetime resolution is 45 fs, almost the same as for the vpK*° decay 
mode. Since the predicted branching ratio for this mode is 0.3%, we obtain 19200 
events in one year of running, with a signal to background of 3:1. Fig. 14 shows 
the x s reach obtainable for a 5cr discrimination between the favorite solution and 
the next best solution, for both decay modes. The background is assumed to be 
20% and the flavor mistag fraction is taken as 25%. The tagging efficiency is taken 
as 10%. The absicca gives the number of years of running, where one year is 10 7 
seconds. The pixel system with the 12 mm gap is called the EOI detector here. 

The other detector configuration that we simulated has the pixels configured 
around the beam leaving a 12 mm x 12 mm square hole. This detector has bet- 
ter efficiency and time resolution (see section V) and now has become the BTeV 
baseline. 

The x s reach is excellent and extends over the entire predicted Standard Model 
range. 



F Measurement of 7 

The angle 7 could in principle be measured using a CP eigenstate of B s decay that 
was dominated by the b — > u transition. One such decay that has been suggested 
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FIGURE 13. The observed decay time distributions for B s — > tpK *° generated with x s = 40. 
Unmixed decays are shown in (a), mixed in (b). Background and mistagging have been included. 
In (c) the results of a likelihood fit to the time distributions are shown. The dashed line shows a 
5cr difference from the best solution. 



is B s — > p°K s . However, there are the same "Penguin pollution" problems as in 
B° — > 7T+7T - , but they are more difficult to resolve in the vector-pseudoscalar final 
state. (Note, the pseudoscalar-pseudoscalar final state here is ir°K s , which does 
not have a measurable decay vertex.) 

Fortunately, there are other ways of measuring 7. CP eigenstates are not used, 
which introduces discrete ambiguities. However, combining several methods should 
remove these. 

We have studied three methods of measuring 7. The first method uses the decays 
B s — > DfK T where a time-dependent CP violation can result from the interference 
between the direct decay and the mixing induced decay [21]. Fig. 15 shows the two 
direct decay processes for B s . 

Consider the following time-dependent rates for neutral B mesons to non-CP 
eigenstates via two different processes that can be separately measured using flavor 
tagging of the other b: 

T(B S -> /) = |M| 2 e _i {cos 2 (xt/2) + p 2 sin 2 (xt/2) - psin(0 + 5) sin(xt)} 

T(B S -> f) = \M\ 2 e- t {cos 2 (xt/2) + p 2 sm 2 (xt/2) + psin(0 — 5) sin(a;t)} 
T(B S -> /) = |M| V*{p 2 cos 2 (>t/2) + sin 2 (a;t/2) - psin(0 - 5) sin(xt)} 
T{B S -> /) = |M| V*{p 2 cos 2 (a;t/2) + sin 2 (a;t/2) + psin(0 + 5) sin(art)}, 
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where M — ( f\B), p — \ y,J , <p is the weak phase between the two amplitudes 
and 6 is the strong phase between the two amplitudes. The three parameters p, 
sin(0 + 5), sin(0 — 5) can be extracted from a time-dependent study if p = 0(1). 

In the case of B s decays where / = D^K~ and / = D~K + , the weak phase is 
7. The decay modes B s -> DfK, D+ -> <f)n + , <p -> i^+i^-, or D s -> K*°K + , 
were simulated. For the 07r + mode, the combined geometric acceptance and re- 
construction efficiency is 5.2% with S/B=10 [22], and the trigger efficiency is 67%. 
In the K*°K + mode the geometric and reconstruction efficiency is 5.9% and the 
trigger efficiencies and signal to background are same as in the 07r + mode. Using 




the branching fractions predicted by Aleksan [21] and assuming a tagging efficiency 
e = 15% we expect 8900 events in 2xl0 7 s. 

The decay time resolution and the detachment of the decay vertex from the 
primary production vertex are shown in Fig. 16 for the D£ — > 07r + decay mode. 
The distributions for the K*°K + mode are similar. 
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FIGURE 
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16. (left) The generated minus measured proper time distribution for B s 
. (right) The distribution in L/a for this mode. 



Using the measured values of S/B and time resolution, a mistag rate of 25%, 
and x s =20, a mini-Monte Carlo was used to generate the extracted value of 7 for 
an ensemble of experiments each with 8900 signal events, for various sets of input 
parameters p, sin(7 + 5),sin(7 — 5). A maximum likelihood fit was then used to 
extract fitted values of the parameters. 

Fig. 17 shows the distributions of the parameters with input values p=0.5, sin 7 = 
0.866 and cos 5 = 0.7. Assuming that sin 7 > then sin 7 can be determined up to 
a two-fold ambiguity, hence 7 up to a four-fold ambiguity. 



50 F 



40 



30 
20 



input value 
Y= 80° 

sin(Y) = 0.866 



0.6 




output 

sin(y) = 0.875 
c= 0.067 

fl 8900 tagged events 



1.1 



1.2 



sin(Y) 

FIGURE 17. Results of determining 7 for many different "experiments" with input value of 
7=80° using the decay time distributions for B° — > DfK T . 



Another method for extracting 7 has been proposed by Atwood, Dunietz and 
Soni [9], who refined a suggestion by Gronau and Wyler [8]. A large CP asymmetry 
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can result from the interference of the decays B — > K D°, D° — > / and £? — > 
K~D , D -> /, where / is a doubly Cabibbo suppressed decay of the _D° (for 
example / = K + tt~, Kittt, etc.) Since B~ — > K~D° is color-suppressed and — > 
K~D° is color-allowed, the overall amplitudes for the two decays are expected to 
be approximately equal in magnitude. The weak phase difference between them is 
7. To observe a CP asymmetry there must also be a non-zero strong phase between 
the two amplitudes. It is necessary to measure the branching ratio B(B~ — > K~f) 
for at least 2 different states / in order to determine 7 up to discrete ambiguities. 
We have examined the decay modes B~ — > K~[K + tt~] and B~ — > if - [fC + 37r] . The 
combined geometric acceptance and reconstruction efficiency was found to be 6.6% 
for the Kit mode and 5.5% for K3ir with a signal to background of about 1:1. The 
trigger efficiency is approximately 70% for both modes. The expected number of 
B ± events in 10 7 s is 2400 in the Kn mode and 4200 in the K3n mode. With this 
number of events we expect to be able to measure 7 (up to discrete ambiguities) 
with a statistical error of about ±8° in one year of running at £ = 2 x 10 32 cm" 2 s~ 1 . 
The overall sensitivity depends on the actual values of 7 and the strong phases. 

The next method, described by Gronau and Rosner [23] and Fleischer and Mannel 
[24], uses B° — > K + n~ and B + K°tt + decays. It is particularly promising as it 
may complement other methods by excluding some of the region around 7 = ir/2. 
We expect to reconstruct 3600 B ± -> K a K± with S/B=0.5 and 29000 B°/B° -> 
K ± ir T with S/B=3. Gronau and Rosner estimate a measurement of 7 to 10° with 
2400 events in each channel [25] , however there has been much theoretical discussion 
about the effects of isospin conservation and rescattering which casts doubt on this 
method [26] [27] [28] [29]. There have been several suggestions, however, on how 
to measure these effects [30], and this method may turn out to be useful. 

VII DECAY TIME RESOLUTION 

In all of these studies we have assumed that we would have 9 /j,m spatial resolution 
in each track hit in the pixel plane. We now address the question of whether or not 
this is reasonable. 

The parameters affecting the pixel resolution include the size of the pixel, the use 
of binary (one bit) versus analog (4 bits) information, threshold and gain variations, 
and the use of electrons or holes as charge carriers, since the drift velocity for 
electrons is three times that for holes. For different incident angles of tracks on the 
pixels the charge sharing is affected by the magnetic field. In BTeV we use a 1.5 T 
dipole field. 

In Fig. 18(a) we show the track angle distribution for two B final states, the two- 
body state 7r + 7r~ and the four-body state n + K + K~ K~ . In both cases the angular 
distribution peaks at small angles, about 50 mr and then falls slowly towards larger 
angles. The spatial resolution has been simulated as a function of angle for various 
pixel sizes. The baseline size is 50 /im x 300 pm. The results of the simulation 
for this size are shown in Fig. 18(b). The magnetic field is in the y direction in 
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this case. So the tracks hitting the x layers are bent. The resolution using binary 
readout is about 10 fim, while it is about 5 /im using 4-bit analog. Note that the 
poorer resolution peak near zero degrees in the non-bend plane does exist in the 
bend plane, but it is shifted toward negative incident angles. Similar results are 
obtained with holes as charge carriers. 
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FIGURE 18. (left) The distribution of track angles with respect to the pixel planes for tracks 
with at least 9 hits for B° — > ir + ir~ and B° — > tt + K~ K + K~ decays, (right) The spatial resolution 
for the X (bend plane) and Y directions, separately, as a function of projected angle, 8^^), for 
binary and 4-bit analog readout using a pixel size of 50x300 /im 2 . 

The resolution in proper time is affected by several factors. One is the inherent 
pixel resolution, as discussed above. Others include the amount of material in the 
pixel system, and the distance the pixel detector is placed from the beam line. 
We take the latter as 6 mm. This distance is limited by the maximum amount 
of radiation damage we are willing to sustain. (The system is retractable during 
machine injection.) In Fig. 19 we show the proper time resolution achievable on 
the if)K*° decay of the B s for several different detector geometries as a function 
of the spatial resolution. The circles represent a geometry with a 12 mm gap and 
equivalent silicon thickness of 600, 500, and 600 /im for the three layers. These 
include the 300 /iin of silicon for each layer, a radio-frequency shield of 100 /zm of 
Al and material for electronics and cooling. 

The simulations presented here have assumed a 9 /im spatial resolution, even 
though we believe that 5 jum is possible. The equivalent silicon thickness of a three 
plane station is taken at 1700 /im. A detector of twice the material thickness and 5 
/iin resolution would have the same time resolution as the one we have been using. 
This points out the need to minimize the material, a well known lesson. 
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FIGURE 19. The calculated proper time resolutions for the decay mode B s — > ipK*°, 
K*° — ► K~ir + as a function of spatial resolution in the pixel system for different detector geome- 
tries and thicknesses. The two geometries considered are a gap of 12 mm between two halves of 
the pixel system and square hole of 12mm x 12 mm. The equivalent silicon thicknesses for the 
three planes are listed. 



VIII COMPARISONS WITH OTHER EXPERIMENTS 

A Comparisons with e + e 5-factories 

Most of what is known about b decays has been learned at e + e~ machines [31]. 
Machines operating at the T(4f>) found the first fully reconstructed B mesons 
(CLEO), B°-B° mixing (ARGUS), the first signal for the b -> u transition (CLEO), 
and Penguin decays (CLEO). Lifetimes of b hadrons were first measured by exper- 
iments at PEP, slightly later at PETRA, and extended and improved by LEP [31]. 

The success of the T(4S) machines has led to the construction at KEK and SLAC 
of two new T(AS) machines with luminosity goals in excess of 3 x 10 33 cm~ 2 s . 
These machines will have asymmetric beam energies so they can measure time 
dependent CP violation. They will join an upgraded CESR machine at Cornell 
with symmetric beam energies. These machines will investigate only B° and B ± 
decays, they will not investigate B s , B c or A& decays. 

Table 4 shows a comparison between BTeV and an asymmetric e + e~ machine 
for measuring the CP violating asymmetry in the decay mode B° — ► 7r + 7r _ . It is 
clear that the large hadronic b production cross section can overwhelm the much 
smaller e + e~ rate. 
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TABLE 4. Number of tagged B° -> tt+tt" (6=0.75 x 1CT 5 ) 





^(cm-^- 1 ) 




# B°/10 7 s 


efficiency 


eD 2 


# tagged 


e + e 


3 x 1CF 


1 nb 


3.0 x 10 y 


0.4 


0.4 


46 


BTeVt 


2 x 10 32 


100/zb 


1.1 x 10 11 


0.06 


0.1 


6400 


BTeV* 


2 x 10 32 


100/xb 


2.1 x 10 11 


0.06 


0.1 


12800 



fThis is for gap detector, expect increase with square hole 
^Luminosity leveled, use 2 x 10 7 s/year (gap detector) 



B Comparisons with Tevatron Central Detectors 

Both CDF and DO have measured the b production cross section [32] and CDF 
has contributed to our knowledge of b decay mostly by its measurements of the 
lifetime of ^-flavored hadrons [33], which are competitive with those of LEP [34] 
and recently through its discovery of the B c meson [35]. These detectors were 
designed for physics discoveries at large transverse momentum. It is remarkable 
that they have been able to accomplish so much in b physics. 

However, these detectors are very far from optimal for b physics. BTeV has been 
designed with b physics as its primary goal. To have an efficient trigger based on 
separation of b decays from the primary, BTeV uses the large 77 region where the 
&'s are boosted. The detached vertex trigger allows collection of interesting purely 
hadronic final states such as 7r + 7r~, D+tt~ and D^K~ . It also allows us to collect 
enough charm to investigate mixing and CP violation. 

The use of the forward geometry also allows excellent charged hadron identifica- 
tion with a gaseous RICH detector. This is crucial for many physics issues such as 
separating Kn from nit, D s n from D S K, kaon flavor tagging etc... 



C Comparison with LHC-B 

LHC-B is an experiment proposed for the LHC with almost the same physics 
goals as BTeV [36]. LHC-B has two advantages: the b cross section is five times 
higher than at the Tevatron while the total cross section is only 1.6 times as large, 
and the mean number of interactions per crossing is three times lower, because the 
LHC has bunches every 25 ns, while the Tevatron bunches come every 132 ns. 

There are, however, many advantages which accrue to BTeV. Let us first con- 
sider the machine specific ones. The 132 ns bunch spacing at the Tevatron makes 
first level detached vertex triggering easier. It is difficult for the vertex detector 
electronics in LHC-B to settle in 25 ns. The seven times larger energy at the LHC 
results in a larger track multiplicity per collision which causes trigger and track- 
ing problems and a larger range of track momenta that need to be analyzed. The 
interaction region at the LHC is relatively short, a— 5 cm, compared with the 30 
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cm long region at Fermilab. This somewhat compensates for the larger number of 
interactions per crossing, since the interactions are well separated. 

There are detector specific advantages for BTeV as well. BTeV is a two-arm 
spectrometer, resulting a factor of two advantage. BTeV has the vertex detector 
in the magnetic field which allows the rejection of high multiple scattering (low 
momentum) tracks in the trigger. Furthermore, BTeV is designed around a pixel 
vertex detector while LHC-B has a silicon strip detector. BTeV can put the detector 
closer to the beam (6 mm versus 1 cm), and has a much more robust tracking 
system which can trigger on detached verticies in the first trigger level, while LHC- 
B triggers on tracks of moderate transverse momentum in their first trigger level. 

We feel that we have more than compensated for LHC-B's initial advantages. 

IX CONCLUSIONS 

Hadron colliders have large b and c cross sections allowing the opportunity for 
precision measurements of CP violation and B s mixing. In our view this requires 
high density tracking and triggering information that can be provided by a state 
of the art pixel system. BTeV has been designed to fit in the new CO interaction 
region at the Tevatron and incorporates a pixel vertex detector, downstream track- 
ing, charged paricle identification, lepton identification and photon detection. The 
vertex detector enables Level I vertex triggering and excellent time resolution on 
heavy hadron decays [37]. 

A summary of the physics reach is shown in Table 5. Those simulations that 
have been upgraded by using the square hole detector are so indicated. 



TABLE 5. BTeV Physics Reach 



Measurement 


Accuracy in 10 7 s 




C = 2 x 10 32 , C leveled 


x s (square hole) 


up to 80 & beyond 


Acp(B" -> 7T+7T-) (gap) 


±0.013 


7 using D S K~ (square hole) 


±8 ot 


7 using D a K~ (square hole) 


±8°* 


B(B~ - K-n+ir) (gap) 


4cr at B of 5.4 x 10~ 8 


sin(2/3) using B° — ► ipK s (square hole) 


±0.013 



f Assumes p=0.7, cos 5=0.7, sin 7=0.5, x s =20 
| For most values of strong phases and 7 



BTeV is an officially recognized R&D project at Fermilab. Development has 
started on the pixel, trigger, RICH, muon, forward tracking and electromagnetic 
calorimeter systems. More information on BTeV can be found on the world-wide- 
web [37]. 
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